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Abstract 



When a small number of individuals of organism of single species is confined in a closed 
space with limited amount of indispensable resources, their breading may start initially 
under suitable conditions, and after peaking, the population should go extinct as the re- 
sources are exhausted. Starting with the logistic equation and assuming that the carrying 
capacity of the environment is a function of the amount of resources, a mathematical 
model describing such pattern of population change is obtained. An application of this 
model to typical population records, that of deer herds by Scheffer (1951) and O'Roke 
and Hamerstrome (1948), yields estimations of the initial amount of indispensable food 
and its availability or nutritional efficiency which were previously unspecified. 
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INTRODUCTION 

The logistic or the Lotka-Volterra model has long been a mathematical frame to 
study population dynamics tending to stationary or oscillating equilibrium due to intra- 
or interspecific interactions (e.g., Pielou, 1974; Begon et al, 1996; Borrelli et al, 1996; 
Glesson and Wilson, 1986; Reed et al, 1996). Also, there is another pattern of population 
change which is singly-peaked. Typical one may be the population change of deer herds 
observed by Scheffer (1951). It was reported that the deer were freed in closed spaces at 
some definite time and the populations increased first nearly exponentially to reach a peak 
and then decreased or went extinct finally. The change was considered to be fluctuation 
or over-abundance from the sigmoidal pattern and ascribed to changes of reproduction 
rate and/or mortality due to unspecified reasons. But, such patterns should be generally 
observable if living organisms are confined in a closed space with constant amount of 
growth resources which are actually not reproducible although initially given. Effects of 
food availability or resource limitation on population dynamics are one of recent concerns 
(e.g., Ogushi and Sawada, 1985; Edgar and Aoki, 1993). To our knowledge, however, 
rather few mathematical models have been studied to analyse such patterns of population 
change and the carrying capacity for population has been traditionally assumed to be a 
constant characterizing its environment. In this report, we propose a new mathematical 
model to interpret such a pattern of population change by introducing a new assumption 
that the carrying capacity is a function of the amount of resources. After formulation and 
its application to the deer herd population, we discuss several characters of our model as 
compared with the existing models. 

MATHEMATICAL MODEL 

We start with the logistic equation for a single species of organism living in some 
limited space, 
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where N is the population size of the organism, r the potential net reproduction rate and 
K the carrying capacity of the population. Now we assume that the carrying capacity 
depends on the amount of indispensable resources in the space for the organisms and the 
resources are consumed by the organisms after they begin to live. Under such situation, 
we may suppose that the carrying capacity is a function of the amount (X) of the resource, 
K = f{X). Then, we have 

1 dN ^ N 

Kl-TT^v) ■ (2) 
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We may further assume that the decreasing rate of X is proportional to the population size 
and the reproduction rate of the resources is negligible compared with the consumption 
rate, i.e., 

dX 

- = -aN (3) 

where a(> 0) is the consumption rate of the resources per individual of the organism per 
unit time. From Eqs.(Q) and (0), we have 

N . s r r x dX 

IxoJiX) 

where No, Xq and to are the initial values of N, X and t, respectively. There may be 
various choices for f(X) as an integrable function which represents a possible resource 
dependence of carrying capacity. We choose here the simplest one, a linear function 
f(X) = bX, with the proportional constant b(> 0), which we may call the nutritional 
efficiency Then we have 

N(t) = No&] r ' ab exp{rt) , (5) 

with t = 0. Equation fl5|) predicts that the amount of the resources per individual, X/N, 
in the case of a = r/b, decreases exponentially with time from the initial value Xq/Nq. 
Solving the simultaneous Eqs.(U) and @, we obtain the following solutions: For the case 
a = r/b, 

N(t) = iVoexp[rt+-^{l-exp(rt)}] , (6) 

r Xq 



and for a ^ r/b, 

N(t) = N [l + (--l)^{l-exp(rt)}Y^ b -^exp(rt) . (7) 
r o A 

The N(t) curve given by Eq.(^) or Eq.(^) has a single peak for a limited range or combi- 
nations of parameters a, b, r, Xq and Nq. The range giving the single peak is determined 
from the extreme condition of N(t). The solution (J^j in the case of a = r/b has the peak 
if rX /aN > 1. We note rX /aN Q = bX /N in this case. The maximum of N is given 
by 

AT rX ° ( aN ° 1\ (Q\ 

N m = exp(— — -1) , (8) 

a rX 

at the time t m = (1/r) ln(rX /aN ). In the case of a ^ r/b, the peak exists again when 
bX /N > 1. The maximum is 

1 J\f T*K T 

iVm . JVo[ _ {r+(a6 _ r) _i }r /(^» ( _^ + 1 __ ) (9) 

with t m = (1/r) \n(rXo/aNo + 1 — r/ab). We show the range where the single peak exists 
on the {^,b) plane in Fig. 1. It should be noted that our model is soluble exactly. We 
also note that it has the scale invariance under the change of (a, 1/6, Xq) into (Aa, X/b, 
XX ) with an arbitrary constant A and the units of X defines the units of a and b. 

APPLICATION TO THE DEER POPULATIONS 

What can be analysed by the present model? To show this, we apply it to the 
population changes of reindeer on St. Paul Island (SPI) from 1911 to 1950 and on St. 
George Island (SGI) from 1911 to 1949 (Scheffer, 1951). The population data have been 
well-known to be of ideal observation in out door laboratory where the animals lived 
under small hunting pressure and were free of predator attack for the 40 years; the definite 
numbers of the animal were planted in the closed spaces at the definite time, after which 
the population showed singly-peaked changes. The accuracy of the numbers was estimated 
to be about 10 %. We also apply the model to the population change of white-tailed deer 



at the George Reserve of the University of Michigan (GRM) which showed a similar trend 
from 1928 to 1947 (O'Roke and Hamerstrom, 1948). 

For the application, we need to fix one of three parameters, a, b and Xq, and need to 
assume the presence of indispensable resources for the animal. We may suppose that it 
was lichen at least for the SPI herd. This is because lichen was considered to be the key 
forage for reindeer, especially in winter (Scheffer, 1951). The grass disappeared on SPI 
40 years after the reindeer introduction, which was regarded as the cause of the reindeer 
extinction. We may apply Eq.([3D here without adding any reproduction term for the 
plant breeding since it was reported that recovery of lichen range may take 15 or 20 years 
there. A caribou is reported to eat 4.5 kg of lichen a day (Bandfield, 1996). We infer that 
real values of the consumption rate of the three deer herds are near to this value since 
they belong to the same family (a Japanese deer is reported to eat 11 kg of grass a day). 
As the choice of the value is not so essential to obtain perspectives to consider the real 
population, we use a = 1.64 tons a year per individual for the three herds commonly. 

The population change (N) of the SPI reindeer from Scheffer's table is shown in 
Fig. 2-A with empty circles. To fit the curve of Eqs. (0), we use the direct search of 
optimization (DSO) for three parameters, r, b and X and obtain r = 0.182 per year, 
b = 0.111 individual per ton and X = 37000 tons. We notice here some deviation of 
the curve from the data points which might be caused by changes of hunting effects or 
weather. We can not clarify the reason at present, however. After the similar application 
of DSO to the population on SGI and that in GRM, the optimized curves are compared in 
Figs. 2-B and C with the observed data. All parameters thus obtained are summarized in 
TAB. I together with the areas of three habitats and the respective initial and maximum 
population sizes. 

Now we explain some characters of the population processes referring the figures and 
the table. The most significant result in the table is that the initial stock X on SPI is 
more than 8 times larger than on SGI although the land areas are almost same. In the 
present model, the deviation of Xq is proportional to that of a due to the scale invariance 
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for parameters mentioned above. However, this difference in X 's is much more than 
one that can be caused by probable difference in a's. Rather, this may correspond to 
about ten times larger N m observed on SPI than that on SGI and suggests that SPI was 
much more fertile than SGI. Scheffer remarked some environmental differences between 
the two islands. Here we propose that the initial values of the carrying capacity is given 
by Kq = bX Q , of which data are also included in TAB. I. K is free from the effect of 
a-ambiguity. The significant difference between Kq's of SPI and SGI in the table also 
supports above view. We find next that the net reproduction rate r of the SPI herd is 
much smaller than that of SGI which is further smaller than that of GRM. Values of r are 
free from the effects of a. A biological reason may exist for the differences of r, although 
we cannot explain it now. The b value of the SPI herd is about twice of that of SGI (and 
GRM). However, this difference might be caused by any difference in possible a values. 

Further, we find significant differences between population processes on SPI and on 
SGI (and in GRM): The population on SPI increased rather slowly and went extinct 
steeply after the maximum while that on SGI increased fast and decayed slowly. For the 
SPI herd, the ratio of the obtained r to the b value is very near to the a value, meaning 
that the curve fitting for SPI reindeer is attained with Eq.(||) or as the case of a = r/b, as 
far as the a value is acceptable. In contrast with this, some similarities are found in the 
population processes of SGI reindeer and GRM white-tailed deer: The r/b is much less 
than the assumed a for both herds, meaning that the fitting is realized with Eq. (|7p or 
as the case of a < r/b. In spite of the large difference between the areas of two habitats, 
the two magnitudes of b are nearly equal each other and the two X 's are too. Between 
the two habitats, a similarity in ecological characters for deers should have existed. 

Now we discuss relations among observed and calculated population parameters. In- 
spection of the table suggests no definite relation of r to N m , X and the respective 
densities, r is presumably inversely related to iVo- The observed N m may have a linear re- 
lation to Xo which is clearly found in Fig. 3. We have shown a non-linear relation between 
N m and X in Eqs. (H) and (H). First, for the case of ab/r = 1, Eq.® is approximated as 
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N m « rX /ea since aN /rX << 1 within the present range of parameters (e is the base 
of the natural logarithm). Second, for the case of ab/r << 1 (the SGI and GRM cases), we 
rewrite Eq.© as N m = N (r/ab) r ^ ab ~ r ^[l-(l-ab/r)N /bX Y^ ab -^ (rX /aN + l-r/ab). 
If ab/r + bX /N » 1, we have N m » (r/a)(r/ab) r/{ab - r) X . This condition is fulfilled 
in the present ranges of the parameters. We have then a linearly-increasing trend of N m 
with X for both cases. Concerned with the coefficients of the linear increase, we show 
the quantitative estimations of the ratio N m /Xo in TAB. II. iV m (DATA) is the maxi- 
mum N, which was really observed on SPI, SGI, or GRM. A m (LNR) is estimated by 
using N m « rX /ea for SPI case or N m « {r / a){r / ab) r ^ ab -^ X Q for SGI and GRM cases. 
Eq.(||) or Eq.(|9]) gives us the fully theoretical value of N m which is denoted by A m (DSO). 
We may consider that these values of coefficients are almost constant over three herds 
causing the linear relation between N m and Xq. 

The minimum requirement of year-long grazing area of lichen for a reindeer was es- 
timated to be 33 acres on SPI (Scheffer, 1951). This meant that the carrying capacity 
per unit area was 0.030 and the carrying capacity of SPI was 800 individuals (Dasman, 
1964). The peak densities (A m /area, estimable in TAB. I) exceeds 0.03 in two habitats, 
SPI and GRM, which were considered to be fluctuations over the carrying capacity. In 
our model, we postulate that the carrying capacity is not a constant of a land but a 
changeable parameter which depends on environmental conditions, e.g., the quantity of 
indispensable forage for the animal. Referring the Kq, the initial value of the carrying 
capacity defined above, we find N m < K in TAB. I, a reasonable limiting relation of the 
maximum population to the maximum carrying capacity. 

Finally, we compare the present model with the original Lotka-Volterra system (LVS) 
for predator -prey interaction. In fact, at a glance, the deer may be regarded as predator 
and the lichen as prey. The system is given by 

dP 

= - cP + a PS , (10) 
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for the predator population size (P) and 

^ = kP-[3PS , (11) 

for the prey population (S) with the coefficients c, k, a and j3 of the well-known meanings 
(Borrelli et al, 1996). When k = 0, this system becomes that of the ordinary differential 
equations of Kermack-Mckendrick type (KMS) and can reproduce a singly-peaked process 
if the initial value of S is larger than c/a. However the LVS or KMS contains as its 
essence the encounter term which is proportional to PS. This means that encounter 
between two interacting species should take place with a constant probability uniformly 
through-out the space and time (applicability of the mass-action law). Hence the system 
should be applicable to the case of thin populations of prey and predator. Our model has 
no such encounter term (see Eq. (|3|)) to lead such limitation to the population density. 
The estimated values of X or X per unit area may be interpreted to be of thin or dense 
population (or stock) of the lichen (or forage) according to its magnitude. For the predator 
or deer, the present model assumes only intraspecific competition as the original logistic 
does. Hence the estimated values of N m and r may be of dense population. Of course, 
effects of overcrowding can be discussed within LVS by introducing the S 2 and P 2 terms 
to it. However, an addition of new terms with new parameters may make the analysis 
more vague unless the parameters are determined by any other methods. We should also 
note that the unimodal curves can be reproduced by a modified logistic equation with a 
term of integrated toxins for population (Small 1987). However, the model has no explicit 
relationship with the resources for the population. 

CONCLUSIONS 

We have presented a simple mathematical model with which one can analyse singly- 
peaked population processes. Although it is simple, the model provides a good account 
of the deer population dynamics by assuming the resource-dependent carrying capacity 
and by introducing two ecological parameters, the consumption rate of indispensable 
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resources (a) and the nutritional efficiency (b), in addition to such traditional ones as 
the reproduction rate (r) and the initial stock of the indispensable resources X . Here a 
and b can be in principle determined by observation. The model is soluble exactly as the 
original logistic is, providing mathematical benefits. It may be applicable to consumption 
of fertilizer by plant (perfectly zero breeding of prey) and to the case of non-zero breading 
of prey by adding a breding term for it in Eq. ([3D. Further we add that a population can go 
extinct steeply, or even suddenly, from its peak in the model. Breeding and extinction of 
many organisms should depend or should have depended on their indispensable resources 
of finite amount to which processes the present model may be applied. 
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TAB. I. Population data of three deer herds: The habitat area, the initial and max- 
imum population sizes (iVo and N m ) are from references (Scheffer, 1951; O'Roke and 
Hamerstrom, 1948). The nutritional efficiency (6), the initial stock of indispensable food 
(Xq) and the reproduction rate (r) which are defined in text are obtained in this work 
after direct search optimization of the theoretical curve (Eqs.(^|) and (0) in text) to fit the 
population records in the references and shown with significant figures of three digits. The 
consumption rate of indispensable food (a) is fixed to be 1.64 tons per year per individual 
for three herds. Kq(= bXo) is the initial carrying capacity given in text. 



Herd 


Habitat area 


N 


N 


r 


b 




K 




(acre) 






Gr 1 ) 


(indiv./ton) 


(ton) 




St. Paul I. reindeer 


26500 


25 


2046 


0.182 


0.111 


37000 


4090 


St. George I. reindeer 


22400 


15 


222 


0.469 


0.0512 


4460 


229 


George Res. w.t.deer 


1200 


6 


211 


0.740 


0.0561 


4150 


233 



TAB. I. 



TAB. II. The estimations of N m /X of three deer herds: A™, (DATA), A m (LNR), and 
A m (DSO) are defined in text. They are divided by X which takes the value corresponding 
to each herd in TAB. I. 



Herd 


A m (DATA) /Ao 


A m (LNR)/A 


A m (DSO)/A 


St. Paul I. reindeer 


0.0553 


0.0407 


0.0409 


St. George I. reindeer 


0.0498 


0.0352 


0.0356 


George Res. w.t.deer 


0.0508 


0.0417 


0.0419 



TAB. II. 
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Figure Captions 

FIG. 1. The range where the condition giving a peak in N(t) curves is fulfilled: 
b > N /X . Notations are defined in text. 

FIG. 2. Population curves obtained with Eqs.(|) and (0) in text to fit the population 
records of the three deer herds. The nutritional efficiency (b) and the initial stock of 
indispensable food (Xq) and the potential reproduction rate (r) which are defined in 
text are optimized. Plate A, reindeer on St. Paul Island (Scheffer, 1951); Plate B, 
reindeer on St. George Island (ibid); Plate C, white-tailed deer in George Reserve 
Michigan (O'Roke and Hamerstrom, 1948). 

FIG. 3. The maximum size of deer populations observed (N m ) vs the initial amount of 
indispensable food (Xq) estimated in text. 
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